Mechanical and dielectric relaxation spectra in seven highly viscous glass formers 
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Published dielectric and shear data of six molecular glass formers and one polymer are evaluated 
in terms of a spectrum of thermally activated processes, with the same barrier density for the 
retardation spectrum of shear and dielectrics. The viscosity, an independent parameter of the fit, 
seems to be related to the high-barrier cutoff time of the dielectric signal, in accordance with the 
idea of a renewal of the relaxing entities after this critical time. In the five cases where one can fit 
accurately, the temperature dependence of the high-barrier cutoff follows the shoving model. The 
Johari- Goldstein peaks, seen in four of our seven cases, are describable in terms of gaussians in 
the barrier density, superimposed on the high-frequency tail of the a-process. Dielectric and shear 
measurements of the same substance find the same peak positions and widths of these gaussians, 
but in general a different weight. 

PACS numbers: 64.70.Pf, 77.22. Gm 
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I. INTRODUCTION 

The publications of Kia Ngai deal with more sub- 
stances and more measurement techniques than the work 
of any other scientist in the field of undercooled liquids. 
Within the past three decades, whenever a new develop- 
ment appeared, he was the quickest to appreciate it, ana- 
lyze it and bring it to the general attention, thus speeding 
up the progress substantially. 

Many scientists in the field share his conviction that 
the flow process in highly viscous liquids can only be 
understood by combining all possible techniques for its 
studj*i^^^. The present paper evaluates recently 
published 7,8 broadband shear and dielectric relaxation 
data on seven glass formers. The mechanical shear data 
were obtained with a new technique 9 which allows to 
cover a large dynamical range. Samples for dielectric and 
shear measurements were taken from the same charge, 
and the temperature sensors of both measurements were 
calibrated to each other. 

The data show a striking similarity of G"(ui) and e"(u>) 
on the right hand side of the a-peak, a similarity which is 
sometimes perturbed by the secondary Johari-Goldstein 
pea k 10 ! 11 . The similarity suggests a common origin of the 
cv-peak in dielectrics and shear. The fact that the shear 
peak appears at a higher frequency than the dielectric 
peak is explainable in terms of the viscosity, which in a 
compliance treatment^ is a free parameter. 

In order to identify the elementary processes with ther- 
mally activated jumps over an energy barrier V, one can 
use a recent translation^ of the textbook^ retardation 
spectrum L(hiT) into a barrier density function l(V). We 
will argue that the compliance barrier density of the shear 
equals the electric dipole moment barrier density l e (V) 
of the dielectric data. 

The next section (section II) explains and motivates 
this approach in more detail. The results of the data 
treatment are described in section III. They are discussed 
and compared to other approaches in section IV. Section 



V summarizes and concludes the paper. 



II. THE BARRIER DENSITY FUNCTIONS FOR 
SHEAR AND DIELECTRICS 

The choice of a retardation spectrum for the shear is 
motivated by a surprising coincidence, which is more or 
less visible in the data of all seven substances. We show 
the two examples where it is most clearly seen in Fig. 1 
and Fig. 2. 

Fig. 1 compares G"(u>) and —e"(u>) at the same tem- 
perature for the type-A glass former PPE. PPE, 1,3- 
Bis(3-phenoxyphenoxy)benzene, is a diffusion pump oil 
with the commercial name Santovac-5P, a molecule con- 
sisting of five phenyl rings connected by four oxygens to 
form a short chain. In terms of the classification proposed 
by the Bayreuth groupie, it is a type-A glass former, a 
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FIG. 1: Comparison of G"[uf) to a properly scaled ~e"(cu) in 
a double- log scale for PPE at 254 K. 
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glass former which shows no or at least no pronounced 
secondary Johari-Goldstein peak. The (negative) e"(a>)- 
data have been scaled to coincide with the G"(w)-data 
on the right hand side of the peak. One finds good agree- 
ment between G"(ui) and —e"{uS) as soon as the frequency 
is two decades higher than the one of the peak in G"{uj). 

In addition, Fig. 1 shows a change of slope of the a- 
tail, from an a; -1 / 2 - to an w~ ^-behavior. The tendency 
is seen most clearly in the dielectric data, but it is also in 
the shear data; their fit improves markedly if one allows 
for a cj^^-component. This might be the influence of a 
hidden Johari-Goldstein peak, an explanation which has 
been favored for glycerol on the basis of pressure, aging 
and chemical series measurements (for a review, see Ngai 
and Paluchi^). But there is also impressive experimen- 
tal evidence for a limiting w _1 / 3 -behavior in shear com- 
pliance data of type-A molecular glass formers^, which 
show the so-called Andradeil creep, J(t) ~ i 1 / 3 , in the 
short-time limit. Therefore we will fit our data in terms 
of a sum of an w _ ^-term (with /3 as fit parameter) and 
an o; _1 / 3 -term, dominating at high frequency. 

Toward lower frequency, the common slope terminates 
for the shear data already at a higher frequency (a shorter 
time) than for the dielectric data. The natural explana- 
tion for this is that the parallel between dielectrics and 
shear is in fact between shear compliance and dielectric 
susceptibility. In a comparison of these two quantities, 
the shear compliance starts to deviate from the dielectric 
susceptibility as soon as the viscous flow sets in. This 
suggests a treatment of the shear data in terms of a re- 
tardation spectrum, with the viscosity as an independent 
parameter—. 

The good agreement at the right-hand side of the a- 
peak is found to be a general feature of all seven mea- 
sured substances, as long as there is no disturbing in- 
fluence from the secondary Johari-Goldstein pea k 10 ' 11 . 
This is seen in Fig. 2, which shows that the good 
agreement between G" and — e" disappears as the peaks 
merge. The substance, tripropylene glycol (TPG), 
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FIG. 2: Comparison of G"{oj) to a properly scaled —e"(ui) in 
a double-log scale for TPG at 188 and 192 K. 



CQH20O4, is still a molecule and not yet a polymer (the 
whole series from the small molecule propylene glycol to 
long-chain polypropylene glycol is well-investigated by 
dielectric a 18 ' 19 ' 20 ). TPG itself has been studied under 
aging 2 - and under pressure 2 ^,. It consists of three con- 
nected propylene groups, with a large dielectric moment 
and a pronounced Johari-Goldstein peak (which is absent 
or at least much less pronounced in propylene glycol^). 

Fig. 2 shows another tendency which will be evaluated 
quantitatively in this paper, namely a much stronger in- 
crease of the imaginary quantities in the a-peak region 
with temperature than in the Johari-Goldstein peak re- 
gion. 

If one wants to decompose a measured relaxation into 
a spectrum of exponential decays in time, one can choose 
between two equivalent possibilities^, the relaxation 
spectrum in which the elementary exponential relaxators 
add to decrease the modulus or the retardation spectrum 
in which they add to increase a susceptibility. In princi- 
ple, the choice is not crucial, because the two spectra can 
be calculated from each other. Here, we choose the retar- 
dation spectrum, with the viscosity 77 as an independent 
variable. 

For this choice, one has the textbook expressions 12 for 
the real and imaginary parts of the complex frequency- 
dependent shear compliance 



j'M = J 9 + 



L(lnr) 



1 + UJ 2 T 2 



d(lnr) (1) 



and 



J» 



i(lnr) 



^d(lnr) + — , (2) 
-~- ojrj 



1 + UJ 2 T 



where r is the relaxation time and L(lnr) is the weight of 
this relaxation time in the retardation spectrum. J g , the 
glass compliance, is the inverse of the infinite frequency 
modulus Goo. 

In an energy landscape picture^ 3 -, one reckons with 
thermally activated jumps over the energy barrier be- 
tween two neighboring minima. In fact, one very of- 
ten finds a broad secondary relaxation peak (the Johari- 
Goldstein pea k 10 ' 11 ) below the a-peak of the flow pro- 
cess. This peak follows the Arrhenius relation in the 
glass phase, indicating that it stems from local thermally 
activated jumps. For a jump over an energy barrier of 
height V, the Arrhenius relation for the relaxation time 
Ty reads 



t v = r e 



V/k B T 



(3) 



where tq = 10 13 s and T is the temperature. 

For a spectrum of thermally activated jumps, one 
defines^ 3 - the barrier density function l s (V) 



ls(V) 



G c 



k B T 



L(V/k B T + ]nT ). 



(4) 
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The index s stands for the shear. With this definition, the 
complex shear compliance equations ([T]) and ^ trans- 
form into 

,00 1 
J'(cu) = J a + J a l s (V) dV (5) 

Jo ' v 

and 

J"{u) = J g l.(V) T l 2 dV + —, (6) 

Jo 1 + UT) 

The dielectric susceptibility can also be describe d 24 i 25 
in terms of a dielectric barrier density function l e (V) 

e(0)- eoo J 'l + uu 2 T 2 

and 

e(0)-eoo J eV 7 1 + w 2 t 2 v ; 

Here e(0) is the static dielectric susceptibility, is the 
real part of e(u>) in the GHz range (larger than n 2 , the 
square of the refractive index, because of vibrational 
contributions^) . 

The above definitions of eqs. ([S][5]) imply a normaliza- 
tion of both l s (V) and l e (V) with 

roc jO _ j 

/ l s{V )dV=±-±, (9) 
JO Jg 

where J° is the recoverable compliance of the steady- 
state flowi 2 -, and 

/>oo 

/ l e (V)dV = l. (10) 
Jo 

The dielectric a-peak occurs always at a lower fre- 
quency than the shear one and seems to coincide 
with the heat capacity and the structural relaxation 
peaksii 2 -^^. Below, we will adopt the view that the 
left side of the dielectric peak marks the disappearance 
and renewal of the relaxing entities. 

In order to describe this decay, one needs to multiply 
the barrier density of the energy landscape with an ap- 
propriate cutoff function at a cutoff barrier V c . Here, we 
will assume that the relaxing entities decay exponentially 
in time with the critical relaxation time r c . With the Ar- 
rhenius relation r c = tq cxp(V r c /fcs7 1 ), this translates into 
a double-exponential cutoff 

c{V) = cxp(- exp((V - V c )/k B T)). (11) 

Equations © and show that a Johari- Goldstein 
peak in G"(lo) or e"(uj) at the peak frequency to\ cor- 
responds to a peak in l(V) at a peak barrier V\ = 
kBT\n(l/ojiTo). We will see that the Johari-Goldstein 
peaks are reasonably well described by gaussians in l(V). 



To describe both the a-peak and the Johari-Goldstein- 
peak in terms of a barrier density, l s (V) and l e (V) will 
be fitted by the form 

l(V) = {a p ^ v ' kBT + a 1/3 e y / 3kBT + ai e~^ v ~ v ^)c(V). 

(12) 

The first two terms describe the high-frequency tail of 
the a-process, the third term the Johari-Goldstein peak 
(if there is one; in three of our seven examples, it is not 
needed). The dimensionless parameter f3 determines the 
slope w _/3 at the beginning of the a-tail in the double- log 
plot of Fig. 1. 

Instead of using the three prefactors ap, ai/3 and a\ 
as fit parameters, it is better to use the correspond- 
ing weights wp, W1/3 and w\ in the integral over the 
barriers, equs. (© and (fTU|) . A type- A glass former 
without Johari-Goldstein peak with w% — is char- 
acterized by the two dimensionless parameters (3 and 
62 = Wi/^/(wf3 + W1/3), at least as far as the form of 
its spectrum is concerned. j3 and 62 have the advantage 
to be reasonably temperature-independent. 

With this prescription, one can fit the e"(u>) of a type- A 
glass former with two temperature-independent parame- 
ters, (3 and b 2 , and two temperature-dependent parame- 
ters, Ae = e(0) — and V c . Their temperature depen- 
dence is a decrease with increasing temperature, which 




FIG. 3: (a) Data and fit of e"(u>) in a double-log scale for 
TPE between 254 and 274 K (b) the same for G(lo). 
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is well fitted by an appropriate power law 



III. DATA EVALUATION 



A<(T) = A l (T, l )(^ r Y\ (13) 



A. The three type-A glass formers 



where T g is the glass temperature. Similarly, one de- 
scribes the decrease of V c with the exponent 7y arL d the 
one of Goo with ^q. 

The strategy of our evaluation is to fit l(V) to the 
dielectric data, and then use the same spectral form 
to describe the shear. The fit of the shear data re- 
quires three additional temperature-dependent parame- 
ters, J g , Jg and r\. Again, it is worthwhile to look for 
combinations which might turn out to be temperature- 
independent. One of them is the ratio 



/o = 



7° - 7 

d e ,J g 

Jn 



(14) 



which appears in the normalization of the shear spec- 
trum, eq. ([9]). A second interesting possibility is not to 
fit the viscosity 77, but the ratio 



foJgV' 



(15) 



where r c is the Arrhenius relaxation time of the terminal 
barrier V c . As we will show in the discussion, one can 
argue that the ratio fj c should be 2 for a renewal of the 
relaxing entities within the critical time r c . 

In the case of a type-B glass former, Fig. 2 shows that 
one needs another dimensionless parameter, because the 
weight of the Johari-Goldstein peak is different in the 
two quantities. In Fig. 2, the Johari-Goldstein peak is 
more prominent in the shear signal, but this varies from 
substance to substance. 
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FIG. 4: Fit values of Goo in DC704 as a function of tem- 
perature. The continuous line is the temperature dependence 
Vc/Av (with Aw = 0.057 ran 3 ) expected from the shoving 
models. 



Three of our seven substances, TPE, DC704 and PPE, 
happen to have no or at least only a rather weak Johari- 
Goldstein peak. Let us begin with TPE. TPE stands 
for triphenylethylene, C2oH\§, a rather flexible molecule 
with three phenyl rings attached to a central C = C 
double bond. Fig. 3 (a) shows data and fit for e"(ui) 
in a double-log plot, Fig. 3 (b) the ones for G(co). The 
dielectric data in Fig. 3 (a) are well fitted with only the 
first two terms of eq. (|12[) . without any Johari-Goldstein 
peak. (3 and hi turn out to be temperature-independent 
within experimental accuracy. 

One gets a good fit for the shear data in Fig. 3 (b), tak- 
ing over (3, b 2 and the cutoff barrier V c from the fit of the 
dielectric data at the given temperature and fitting Goo, 
/o and fj c . Goo is temperature-dependent, but /o and 
fj c are again temperature-independent within the exper- 
imental accuracy, thus justifying our choice of variables. 
The parameters and their temperature dependence are 
listed in Table I. 

The temperature exponents 7y and 70 of the critical 
barrier V c and the infinite frequency shear modulus Goo 
are the same within their error bars (about 5 % for 7y 
and about 10 % for 70). This shows the validity of the 
shoving models, according to which the energy barrier 
of the a-process should be proportional to the infinite 
frequency shear modulus Goo ■ The shoving model postu- 
lates that the a-process happens when the local energy 
concentration exceeds the product GooA„, where A„ is a 
volume expansion. 

The same results, maybe even a bit clearer because of 
the stronger dielectric signals, are obtained for the two 
other type-A glass formers PPE and DC704. Again, the 
fit parameters are listed in Table I. In particular, the 
cj~V 3 -contribution j s mucn better seen, as illustrated 
in Fig. 1 for PPE. In DC704, again a diffusion pump 
oil (l,3,3,5-tetramethyl-l,l,5,5-tetraphenyl-trisiloxane, a 
rather large molecule) we have the additional advan- 
tage of a large temperature range of the measurement, 
from 209 to 239 K. As in TPE, we find temperature- 
independent parameters 0, 62, /o and fj c . Again, we 
find the shoving mode&. confirmed in both glass form- 
ers. In DC704, one even sees the curvature of both curves 
(see Fig. 4), which justifies our temperature exponent 
Ansatz, eq. (fT3")l . 

Table I comprises the fit parameters for these three 
type-A glass formers. Note that our formalism allows to 
describe both shear and dielectric data over the whole 
temperature range with eleven temperature-independent 
parameters. 
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glass former 


TPE 


DC704 


PPE 


T a (K) 


249 


211 


244 


Ae 


0.0491 


0.257 


2.011 


7e 


1.85 


2.26 


1.90 




0.77 


0.85 


1.04 




0.18 


0.27 


0.215 


W) (eV) 


0.767 


0.639 


0.755 


-yv 


4.3 


4.6 


4.6 


G x (T g ) (GPa) 


1.38 


1.80 


1.27 


7G 


4.7 


4.2 


4.7 


fa 


1.65 


2.38 


2.22 


fjc 


2.5 


2.45 


2.05 



TABLE I: Parameters of the three type-A glass formers. Up- 
per part e(w), lower part G(lo). 



B. The four type-B glass formers 

In the type-B glass formers DHIQ, PB20, Squalane and 
TPG, one needs to fit a Johari-Goldstein peak on top of 
the high-frequency tail of the a-process. This is illus- 
trated in Fig. 5 for our first type-B example, squalane. 

Squalane, C^H^, is a short chain molecule with 24 
carbon atoms in the backbone and 6 attached CH 3 - 
groups, rather polymerlike. It has a strong and well- 
separated Johari-Goldstein peak (see Fig. 5), much bet- 
ter visible in the shear data than in the dielectric data. 
The dielectric dipole moment is very weak. Nevertheless, 
it is possible to fit both sets of data with the same retar- 
dation spectrum, attaching a substantially higher weight 
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FIG. 5: Data and fits of (a) e"(w) (b) G(ui) in squalane. 



to the Johari-Goldstein peak in the shear (see Table II). 

In this substance, it is not possible to fit the shear data 
with a temperature- independent parameter /o; one has 
to postulate a rather strong increase of /o with increasing 
temperature 

f (T) = f (T g )+f^(T-T g ), (16) 

but one can keep the parameter fj C constant (see Table 
II). 

PB20 is a true polymer, relatively short (5000 g/mol), 
composed of 80 % 1,4-polybutadiene monomers and 20 
% 1,2-polybutadiene monomers. The results look very 
similar to those of squalane, and the resulting fit param- 
eters in Table II are in fact close to those of squalane. 
Even more than squalane, it has the polymer feature of 
a relatively slow decrease of e '(u>) at low frequency, ex- 
plainable in terms of chain modes with long relaxation 
times^. This is illustrated in Fig. 6, which shows the 
deviation between fit and data at low frequency. As a 
consequence, the resulting parameters have a larger er- 
ror bar in squalane and polybutadiene than in the two 
molecular substances TPG and DHIQ. In particular, the 
deviations between and 7g do not demonstrate a fail- 
ure of the shoving model. 

TPG is a much more favorable case, with a very large 
dipole moment and no problems at the cutoff barrier. As 
Fig. 7 (a) shows, our spectrum of eq. (fF2|) provides beau- 
tiful fits over a large temperature range. One needs to 
take the temperature dependence of the Johari-Goldstein 
peak position V\ into account. Our fit found 

Vl = 0.295^, (17) 

a bit smaller shift than the one found in aging 
experiments^. 

Fig. 7 (b) shows that the shear data are well described 
in terms of the dielectric retardation spectrum. There is 
a small temperature dependence of /o, but fj C is again a 
temperature-independent constant. The shoving model 
is found to be well fulfilled (see Fig. 8). 
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FIG. 6: Data and fit of e"(u;) in polybutadiene at 180 K. 
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Finally, DHIQ, decahydroisoquinoline, CgHnN, is 
best described as two cyclohexanol rings sharing one 
C — C-bond, one of the two rings having an NH replacing 
a Ci/2-group. In this case, the Johari-Goldstein peak is 
very prominent in the dielectric dat a 27 i 28 , comparable to 
the one in G(ui). The dipole moment is large; both V c 
and Goo can be determined with high accuracy. Again, 
their temperature exponents 7y and agree within the 
error bars (see Table II), in agreement with the shov- 
ing model 2 ^. Since both are exceptionally large (DHIQ 
is very fragile, m=158 in Angell's scheme 2 ^), their good 
agreement provides a strong argument for the validity of 
the model. 

In Table II, the Johari-Goldstein peak is characterized 
by the weight of the peak 

w(T) = aiV^Tr^x) = ai FWHM^f{Tt/A\n2) (18) 

which shows a Boltzmann factor behavior 

w{T) = w{T g ) exp(-E a (l/k B T - l/k B T g )), (19) 

with a formation energy E a which is on the average 2/3 
of the peak position V\. 



IV. DISCUSSION 

The preceding section presented a quantitative descrip- 
tion of the a- and the /3-process in dielectrics and shear 
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glass former 


Squalane 


PB20 


TPG 


DHIQ 


T 3 (K) 


167 


176 


184 


175 


Ae(T 9 ) 


0.0155 


0.132 


23.3 


1.707 




2.1 


0.0 


1.51 


0.0 





0.6 


0.44 


0.85 


0.4 




0.2 


0.2 


0.22 


0.2 


Vc(T g ) (eV) 


0.517 


0.54 


0.63 


0.635 


"IV 


3.2 


3.8 


3.0 


6.4 


G 00 (T 9 ) (GPa) 


1.33 


1.63 


2.69 


3.1 


7G 


2.5 


2.7 


2.8 


6.3 


fo(T g ) 


2.4 


3.25 


6.7 


1.56 


fo (1/K) 


0.4 


0.41 


-0.04 


0.25 




2.7 


2.4 


2.5 


2.0 


Vx (eV) 


0.27 


0.28 


0.32* 


0.32 


FWHM(eV) 


0.135 


0.170 


0.154 


0.16 


W S {Tg) 


0.56 


0.55 


0.15 


0.50 




0.03 


0.23 


0.02 


0.48 


E a (eV) 


0.25 


0.12 


0.19 


0.25 


*average value, see eq. 











TABLE II: Parameters of the four type-B glass formers. Up- 
per part G(u), middle part t(uj), lower part Johari-Goldstein 
peak parameters for both. 



for seven different glass formers, a description which is 
based on the concept of isolated and independent ther- 
mally activated jumps in the energy landscape. The de- 
scription allows for a reasonable fit of the temperature 
dependence in terms of temperature- independent param- 
eters. The number of parameters is not small; one needs 
eleven or twelve parameters for a type-A glass former 
(depending on whether /o is temperature-independent or 
not, see Table I and II) and five additional parameters 
for the description of the f3- or Johari-Goldstein peak (see 
Table II). 

Nevertheless, the exercise is not completely meaning- 
less. One does indeed get meaningful quantitative in- 
formation, which is impossible to obtain otherwise. The 
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FIG. 7: Data and fits of (a) e" (u) (b) G(w) in TPG. 
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FIG. 8: Fit values of Goo in TPG as a function of temperature. 
The continuous line is the temperature dependence V c /Av 
(with Av = 0.037 nm 3 ) expected from the shoving model 2 ^. 
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first and rather important one is the probable equality of 
the retardation spectra of shear and dielectrics (but with 
a different weight of the Johari-Goldstein peak) , an infor- 
mation which one can guess from the raw data (see Figs. 
1 and 2), but which requires a full fit for its quantitative 
check. 

The second and equally important quantitative infor- 
mation concerns the dimensionless ratio fj C between the 
terminal dielectric relaxation time r c and the product 
of the viscosity with the total retardation compliance, 
eq. (fT5|) . The seven fitted values lie between 2 and 2.7 
(average value 2.37). This indicates a general relation 
between the dielectric terminal time and the viscosity. 
Since the dielectric terminal time seems to coincide with 
the structural lifetime^, it is probably also the lifetime of 
the double-well potentials which are responsible for the 
retardation spectrum. 

Question: What do we expect for the ratio fj C if this 
is indeed the case? To answer this question, consider 
a constant applied shear stress. After the time r c , all 
the double-well potentials of the spectrum would have 
reached thermal equilibrium, giving their full contribu- 
tion to the compliance. From this consideration, if we 
renew them at the time r c , we would naively expect them 
to be able to give their contribution again after this time, 
yielding fj c = l. 

But this answer is not correct. To get the correct an- 
swer, one must consider the difference between energy 
and free energy in these double- well potentials. To keep 
the argument simple, let us restrict ourselves to the spe- 
cial case of a symmetric double- well; it applies as well to 
the asymmetric case. 

If the double- well is initially symmetric and if it couples 
to the shear stress a with a coupling constant v (the 
coupling constant has the dimension of a volume), then 
the asymmetry A under the stress is av. One well has 
the energy —av/2, the other has the energy +av/2. 

In thermal equilibrium, the population of the two wells 
is given by their Boltzmann factors. It is easy to calculate 
the energy U of the equilibrated system in the limit of a 
small stress 



U 



(20) 



This is the energy transported to the heat bath in the 
equilibration of the relaxing entity after switching on the 
stress. 

The free energy F is 



F 



(21) 



only half of the energy itself. If one thinks about it, the 
reason is clear: spending the energy, one has spanned an 
cntropic spring by the population difference in the two 
wells. If one removes the stress slowly, one gets half the 
energy back. But if the double-well potential decays, one 
gets nothing back. 



The contribution of the relaxing entity to the compli- 
ance is given by the second derivative of the free energy 
with respect to the stress. But if we now deal with the 
effect of a renewal of the double- well potential on the vis- 
cosity, we have to count the energy. This means we spend 
twice as much energy under a constant stress as the one 
calculated above in our first oversimplified picture. And 
this means the viscosity must be a factor of 2 smaller, 
which implies fj c — 2. This is reasonably close to the 
fitted values in Table I and II. 

A third quantitative conclusion of the present study 
is a surprising agreement with the conclusions of Plazek 
et ali£ from their recoverable shear compliance experi- 
ments. If one takes the parameters of Table I to calculate 
the recoverable compliance, one gets curves which closely 
resemble those reported by them. Obviously, it is exper- 
imentally much easier to detect the Andrade creepii in 
creep experiments than in dynamical ones. If one calcu- 
lates /o from their data, one finds values between 1.5 and 
2.3, similar to those in Table I. 

Here, however, a word of caution is in place. Our data, 
taken as they are, do not imply a limited recoverable 
shear compliance. In fact, they are well fitted by the BEL 
models, which has a divergent recoverable compliance. 
The values in the two tables stem from the assumption 
that the two retardation spectra of dielectrics and shear 
(at least as far as the a-peak is concerned) are the same. 

The same is true for the fourth conclusion, the valid- 
ity of the shoving model^. The fitted Goo-values were 
obtained under the same assumption. 

Finally, the Johari-Goldstein peak increases its height 
with increasing temperature. The increase follows a 
Boltzmann factor, with a formation energy of about two 
thirds of the barrier height at the center of the peak. 



V. SUMMARY AND CONCLUSIONS 

Dielectric and shear relaxation data in seven highly 
viscous liquids, most of them molecular liquids, were 
evaluated in terms of a barrier density of independent 
thermally activated relaxation centers. Three of the sub- 
stances are type-A glass formers without or with only a 
very small Johari-Goldstein peak, four of them show a 
pronounced Johari-Goldstein peak. 

The most important conclusion is the probable equal- 
ity of the dielectric and shear retardation spectra, 
guessed from the raw data and confirmed by a quanti- 
tative fit. The difference in the peak positions is due to 
the influence of the viscosity. The Johari-Goldstein peak 
has different weight in dielectrics and shear. 

The second important conclusion concerns the viscos- 
ity. It seems probable that the viscosity results from 
the constant renewal of the double-well potentials in the 
sample within the terminal dielectric relaxation time. 

Our data support earlier recovery compliance results 
by Plazek et aU^, according to which one has an 
Andrade^ creep J ~ t 1 / 3, at short times in type-A glass 
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formers (glass formers without Johari-Goldstein peak). 
They further support the shoving models, which pos- 
tulates a proportionality between the infinite frequency 
shear modulus and the Arrhcnius barrier of the terminal 
relaxation time. 
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